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Abstract. In this paper, we fully investigate the cosmological effects of the 
moduli-dependent one-loop corrections to the gravitational couplings of the string 
effective action to explain the cosmic acceleration problem in early (and/or late) 
universe. These corrections comprise a Gauss-Bonnet (GB) invariant multiplied 
by universal non-trivial functions of the common modulus a and the dilaton 
4>. The model exhibits several features of cosmological interest, including the 
transition between deceleration and acceleration phases. By considering some 
phcnomenologically motivated ansatzs for one of the scalars and/or the scale 
factor (of the universe), we also construct a number of interesting inflationary 
potentials. In all examples under consideration, we find that the model leads 
only to a standard inflation {«> > — 1) when the numerical coefficient 5 associated 
with modulus-GB coupling is positive, while the model can lead also to a non- 
standard inflation (w < —1), if <5 is negative. In the absence of (or trivial) coupling 
between the GB term and the scalars, there is no crossing between the < — 1 
and ui > —1 phases, while this is possible with non-trivial GB couplings, even 
for constant dilaton phase of the standard picture. Within our model, after a 
sufficient amount of e-folds of expansion, the rolling of both fields </> and a can 
be small. In turn, any possible violation of equivalence principle or deviations 
from the standard general relativity may be small enough to easily satisfy all 
astrophysical and cosmological constraints. 
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1. Introduction 

Dark energy, or cosmological vacuum energy, spread on every point in space exerts 
gravitationally smooth and negative pressure and hence leads to cosmic inflation, a 
phenomenon which is not accounted for by ordinary matter and General Relativity 
(GR). Unfortunately, at the moment the origin of dark energy (and its associated 
cosmic acceleration) is arguably the murkiest question in cosmology - and the one 
that, when answered, may provide deepest insights into the origin (and the fate) of 
our universe, establishing synergies between the microscopic physics and cosmological 
length scales. 

In recent years, the cosmic acceleration problem has been addressed in hundreds of 
papers proposing various kinds of modification of the energy-momentum tensor in the 
'vacuum' - the source of gravitationally repulsive negative pressure. Some examples 
of recent interest include brancworld modifications of Einstein's GR , scalar-tensor 
gravity 0, phantom (ghost) field [2], K-essence ^ or ghost condensate 0, f{R)- 
gravity that adds terms proportional to inverse powers of i? [Hj or i?" {n > 2) terms, 
or both these effects at once [7|, thereby modifying in a very radical manner the 
Einstein's GR itself. All these proposals typically include a more or less disguised 
feature, non-locality, other than classical or quantum instabilities. A recent review by 
Copeland, Sami and Tsujikawa [H| summarizes the prospects and limitations of these 
proposals in an elegant way. 

String theory is much more ambitious and far reaching than any of the above 
ideas, and it may be compatible with a dc Sitter universe (see, e.g. |H1). It has 
been known that a cosmological compactification of classical supergravity - the low 
energy description of string/M theory - on compact hyperbolic spaces ^3 o^' 
Ricci-fiat twisted spaces [l2|, naturally leads to a period of cosmic acceleration. 
Furthermore, fundamental scalar fields abundant in such higher-dimensional gravity 
theories potentially provide a natural source for gravitationally repulsive vacua, as 
well as a mechanism for generating the observed density perturbations In this 

respect, the best framework where one can hope to fully explain the dark energy effect 
is perhaps string/M theory. 

A strong motivation for a cosmology based on string theory ^1], that goes 
beyond the standard cosmology and incorporates Einstein's theory in a more general 
framework, is that it may shed light on the nature of the cosmological singularities 
arising in GR. Moreover, string theory has the prospect of unification of all other 
interactions, see |15) for a review. An interesting question is which version of string 
theory, if any, among the presently known five (which are related to each other by 
various dualities ^HDj best describes the physical world. A successful superstring 
theory should naturally incorporate realistic cosmological solutions. 

String theory, as a theory of quantum gravity in one time and nine space 
dimensions, literally contains (or predicts) an infinite number of scalar fields. Most 
of the scalars in string theory, arising from a compactification to four-dimensions, are 
structure moduli associated with the (six-dimensional) internal geometry of space 
(Calabi-Yau spaces or Calabi-Yau shapes), which do not directly couple with the 
spacetime curvature tensor and hence to Einstein gravity |17) . It is well appreciated 
that, at least, one of the Kahlcr moduli (or the modulus associated with the overall 
size of the internal compactification manifold) naturally couples to (standard) Einstein 
gravity via Riemann curvature invariants |18l 119) . 

According to conventional wisdom superstrings involve a minimum length scale 
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only a few order of magnitude larger than the Planck length, Ip ^ 10^ cm, and, 
in turn, string theory leads to a modification of Einstein's theory only at very short 
distances. This idea is perfectly reasonable when we study small fluctuations in the 
fixed (time-independent) background, but may not be generally true, for the scalars 
such as the dilaton (0) and the modulus field (cr) and their couplings to various 
curvature terms, which arise universally in any four-dimensional heterotic superstring 
model, may be important in any cosmological solution of the string effective action 
(see, for example [201 1^)- In this paper, we therefore investigate as to whether the 
evolutions of (f> and cr, and their couplings to a Gauss-Bonnet term, can explain the 
cosmic acceleration problem at early (and/or late) universe. 

This paper is organized as follows. In section 2 we briefly discuss the limitations 
of a standard cosmological model where Einstein's gravity is coupled to a scalar field 
with a self-interaction potential alone. In section 3 we first motivate ourselves to 
a more general action which consists of, in addition to a standard scalar potential, 
scalars dependent loop corrections to the gravitational couplings of the string effective 
action and write the full set of field equations for the metric, the modulus and the 
dilaton. In section 4 we first write down the equations of motion as a monotonous 
system of first order differential equations, and present various analytic solutions by 
allowing at most two variables of the model to take some specific forms. In section 
5, we do a similar analysis but expressing the field equations in an alternative form, 
which we find more appropriate for implementing our knowledge about the coupling 
functions and the behaviour of the accelerating solutions at early or late-times. At the 
end of both sections 4 and 5 we provide an overview of the different solutions being 
discussed there. In section 6, we will have a cursory look in the system where gravity 
is non-minimally coupled to mater and radiation. Our conclusions are summarized in 
section 7. 



2. The standard scenario 

In a standard cosmological scenario, one assumes that for a given scalar field a with 
a self- interaction potential V{a)^ the Einstein gravity is described by the action: 

1 

2~2" 2' 



S = 



(1) 



where k is the inverse Planck mass nip^ = {SttGnY^'^ , 7 is the coupling constant 
and a is the classical scalar field whose stress-energy tensor is set by the choice of 
the potential. What one normally wants is a well motivated potential for the desired 
evolution. This often refers to a cosmological potential that allows at least an epoch 
of cosmic infiation or acceleration/deceleration regime. 

A real motivation for an effective action of the form arises from the fact that 
inflation with the dynamics of a self-interaction potential, V{a) > 0, provides a natural 
mechanism for generating the observed density perturbations |13j . To this end, one 
defines the spacetime metric in standard Friedmann-Robertson- Walker (FRW) form: 

3 

ds^ = ^dt^ + a{tf"^{dx'f, (2) 

i=l 

where a{t) is the scale factor of the universe. The Hubble parameter that measures 
the expansion rate of the universe is given hy H = d/a, where the dot denotes a 
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derivative with respect to proper time t. One also assumes that the scalar field a 
obeys an equation of state (EoS), with an equation of state parameter given by 

where q = —ad/a^ is the deceleration parameter. This definition oiw does not depend 
on the number of (scalar) fields in the model. The 'equation of state' parameter w 
is essentially a measure of how squishy the substance is. It is when w < —1/3 a 
gravitationally repulsive pressure exerted by one or more scalar fields can drive the 
universe into an accelerating phase {q < 0). In particular, the cosmological constant 
A as a source of dark energy or vacuum energy means that wa = pa/ PA = ^ 1- 

In the last two decades, many authors have studied cosmology based on (Q) 
by considering various ad hoc choice of the scalar potential V{a) (the recent 
reviews [55] provide an exhaustive list of references). Written in terms of the following 
dimensionless variables: 

X^.^l(^)\ y^.^^, e^A, (4) 

the field equations, obtained by varying the action satisfy the following simple 
relationships 

y = 3 + e, X^-e. (5) 

The time-evolution of the scalar field is given by a{t) = K-i/di[-2ij/7]i/2. Clearly, 
the model lacks its predictability, since e (and hence 3^) is arbitrary. This arbitrariness 
in the potential has motivated many physicists/cosmologists to consider models with 
some specific choices of the potential, e.g. quadratic/chaotic potential V{a) = 
Vq + ^m^tr^, power-law potential V{a) = ^ mj^^'^Xn cr", inverse power-law potential 
V{a) = A4(cro/cr)" {n > 2), the axion-potential Via) = A^{C ± cos((t/(To)), etc. This 
approach to the problem, starting from Linde's model of chaotic inflation j24l . may 
be motivated from different particle physics models or effective field theory. Without 
a complete theory of fundamental gravity, the validity of these potentials (to explain 
the dark energy problem) is not clear. The choice of possibilities for the potential also 
mimics the fact that the exact nature and origin of an inflationary potential (or dark 
energy) has not been convincingly explained yet. 

A common (if not conservative) argument sometimes expressed in the literature 
is the following. To explain a particular value of one new parameter (acceleration, 
i.e., the second derivative of the scale factor, a), it is enough to introduce another 
new parameter in the set up of the model. To this aim, one often fixes the time 
dependence of the scale factor, by making a particular ansdtz for it, which then fixes 
the parameter e, and then finally find a potential, using Q-®, which supports this 
particular a{t). For instance, for the power-law expansion, a{t) cx t", such a potential 

is given by F = Vb e^V^/" (c'/mpi)^ w^is fashion one may find a good approximation 
at each energy scale, but the corresponding solution may have little relevance when 
one attempts to study a wider range of scales, including an early inflationary epoch. 

Action is remarkably simple but not sufficiently general as it excludes, 
for example, the coupling between the field a (a run away Kahler modulus) and 
the Riemann curvature invariants, as predicted correctly by some version of string 
theory Also the best motivated scalar-tensor theories which respect most of 

GR's symmetries naturally include the GB type curvature invariants multiplied by 
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field dependent couplings. That is to say, the most general scalar-tensor theories may 
involve one more parameter, the GB integrand multiplied by a non-trivial function of 
a or (p. One may wonder whether such extra degree of freedom added to the standard 
model ruins the theory, leading to a plethora of both qualitative and analytic behaviour 
of cosmological solutions. This is generally not the case; the knowledge about the a 
(and/or 0)-dependent couplings may even help to construct a realistic cosmological 
model with time- varying e (which otherwise remains arbitrary) without destroying the 
basic characteristics of the standard model 



3. Effective action and equations of motion 



String theory gives rise to two kinds of modifications of the Einstein's theory. The 
first one is associated with the contribution of an infinite tower of massive string 
modes and leads to a'-corrections, and the second one is due to quantum loop effects 
(string coupling expansion) |18j . Therefore, if we wish a cosmological model to inherit 
maximum properties of the low energy string effective action and its associated early 
universe cosmology, then it may be necessary to include an additional scalar, namely, 
the dilaton field 0, which plays the role of string loop expansion parameter. To this 
aim, the effective action for the system may be taken to be ^| EO] 



±.R _ I(Va)2 - V{a) - ^-{V4>f + \[m - 6a^)]nl_ 



(6) 



where 7^^5 = - 4i?^^i?^'' -I- Rf.^paR'"''"' is the Gauss-Bonnet (GB) integrand. 
It is assumed that there is no potential term corresponding to the field (f>. The GB 
term TZq^ has been multiplied by non-trivial functions of (p and a, otherwise it would 
remain a topological invariant in four dimensions. One would expect TZq^ to be small 
in the present epoch and that mass scales associated with the dilaton and the modulus 
to be smaller than mpi, so that they would not be dominating the dynamics. Only 
when these two conditions are met, can the higher-derivative corrections to the action 
be neglected and does string theory go over to an effective quantum field theory )14j . 
Furthermore, we will consider for most of this paper that 7 > 0, C > 0, so that both 
(7 and (/> are canonical (real) scalar fields. 

It is worth noting that a GB integrand multiplied by a function of (j) is already 
present at the string tree- level in the next to leading-order a' expansion ^Sj , and a non- 
trivial ^((t) exists at the one- loop order ^H]- One also notes that the (semi)classical 
tests of Newtonian gravity typically deal with small fluctuations in the fixed (time- 
independent) background and thus are unaffected by the GB type modification of 
Einstein's theory. Another beauty of the action © is that the equations of motion for 
the modulus-dilaton-graviton system form a set differential equations with no more 
than second derivatives of the fields or the metric, which can be solved exactly with 
some specific choice of the (coupling) parameters, or a prior knowledge about the 
evolution of one of the scalars, and/or the scale factor, a{t). 

As compared to the four-dimensional perturbative string effective action obtained 
from a compactification of lOd heterotic superstring on a symmetric orbifold |19| . 
action © lacks the terms proportional to RR = rj^^'^'^^Rfj.u'^'^RpXaT (where rj'^'^P^ = 
e^^P^ I ^—g and a totally anti-symmetric tensor) which however give a trivial 

contribution in the background Q. The numerical coefficient b typically depends on 
the masslcss spectrum of every particular model; in jl9) . it was found to be associated 
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with the four-dimcnsional trace anomaly of the N ~ 2 supcrsymmctry of the theory. 
Here we absorb this coefRcient into the couphng function £,{(j). 

An often asked question is the following: what are the choices for V{a), ^(cr) 
and X{4>))7 And what could motivate these choices? One can take a different 
point of view: it may be more productive to know what form of the potential (or 
the scalar-GB couplings) offers the best resolution to the dilemma posed by the 
current cosmic acceleration or the dark energy problem. In some perturbative string 
amplitude computations, V{<7) appeared to be absent implying that it is a 

phenomenologically- motivated term. Whatever the origin of the potential, V{a), in 
((HJ may be, we believe that the cr-dependence of the potential is quite universal. This 
was also the case revealed recently from type IIB string theory compactified on a 
(deformed) conifold [HI; this study also includes some supersymmetry breaking non- 
perturbative effects, mainly those coming from different branes and fluxes present in 
the extra dimensions. Typically, the moduli potentials are generated by fluxes, so 
V{a) in (|SJ) may take into account a contribution from fluxes. 

Recently, examined a cosmological solution for the system (jHj) with the 
simplest choice of the scalar potential and the modulus-GB coupling, namely V{(t) = 
Voe~'^'^^'^° , ^{a) = ^oe'^'^/'^°, deleting the ^-dependent terms. Even this simple 
system led to plethora of qualitative behaviors of cosmological solutions, among 
which one may accommodate lots of data, including the dark energy's 'equation of 
state' parameter, w. However, a solution of the type a{t) oc t" with a > 0, or 
a{t) oc {too — t)°' with a < 0, i.e. a power-law expansion for which e = const, does 
not support the transition from acceleration to deceleration as well as the crossing 
between quintessence {w > —1) to phantom {w < —1) phases. A non-constant e 
is generally required for achieving a natural exit from inflation in the early universe 
cosmology as well as for explaining a transition from deceleration to acceleration in the 
expansion of the (late) universe, which has been viewed as a recent phenomenon |26| . 
Thus we shall normally consider a situation in which e (and hence the EoS parameter, 
w = — 2e/3 — 1) dynamically changes as it happens in inflationary cosmology. 

A study has been recently made in |27] by considering a variant of the action © , 
which includes additional (/(-dependent higher curvature terms and also a non-trivial 
coupling between cj) and the Einstein-Hilbert term, i.e. in the presence of a single 
scalar field (a dilaton or compactification modulus). The method authors used there 
to find the solutions was mainly to identify analytically some asymptotic solutions and 
then try to join them numerically. A similar approach was taken previously in the 
paper [20], for the action © with V{cr) = ([2111201 provided further generahzations 
with a non-zero spatial curvature). In this work, we examine the system analytically, 
by writing the field equations as a monotonous system of first (and second) order 
differential equations. In the background |(2J), the equations of motion derived from 
the action (O are given by (see, e.g., [201130!) 



= -^H^ + ^&' + U^ + V{a)-3H'f, (7) 



0^,(. + 3i/.) + ^^ + i^7^^„ (9) 

o = C(0 + 3i^0)-i^7^^^, (10) 



On compatibility of string effective action with an accelerating universe 



7 



where / = \{<t>) — C(o'). Note that = const is a solution of the last equation above 
only if A((/)) ~ const, implying a trivial coupling between (j) and the GB term. Because 
of the Bianchi identity, one of the equations is redundant and may be discarded, as 
it will be trivially satisfied by a solution obtained by solving three equations only: 
equations Q and ((Sj) arc not functionally independent. More precisely, the linear 
combination iJ x (7) + iJ x (8) — ^ a x (9) — i x (10) yields the time derivative of (7). 
One may drop, for instance, (8), so that (7), (9) and (10) form a set of independent 
equations of motion. 

In |23| . we analysed the model by treating as a constant, which then leaves 
only cr as a dynamical field. This analysis was partly motivated from an earlier 
observation by Antoniadis et al. [201 j in the V{a) = case, that the behaviour of 
solutions depends crucially on the form of cr-dependent string one-loop corrections, 
while the ^-dependent contribution is negligible and may be ignored. This has got 
another clear motivation: in string theory context, fixing of the volume (or Kahler) 
moduli, is only an approximate idea (see, e.g.pi^). 

We are not interested in the situation where the contributions coming from 
the GB term dominates the other terms in the effective action, which is consistent 
with our approximation of neglecting higher derivative terms, like (Vct)'' and also 
higher curvature terms. For all solutions to be discussed below, the squared of the 
Hubble parameter, H"^ , will decay faster than the couplings A(^) and ^(ct) grow with 
time. Thus the moduli/dilaton dependent two (and higher) loop corrections to higher 
derivative terms or terms higher than quadratic in the Riemann tensor, which give 
at least sixth powers of i?, are small and may be neglected. Nevertheless, it is 
interesting to note that a four-dimensional string effective action without a Gauss- 
Bonnet term, but including terms up to quartic in cr, (Vcr)^, give rise to solutions 
supporting inflation and/or dark energy universe. There also exist models where an 
accelerated expansion can be realized without a field potential; some examples are the 
K-inflation and dilatonic ghost condensate |32| . 

We study here a model without matter couplings; only the matters that are 
present in a fundamental theory such as a string/M theory arise from the decay of the 
fields, like a and </>. In the last section, however, we will have a cursory look in the 
system where 4> is coupled (non-)minimally to ordinary matter and radiation. 

4. First-order system 

Let us define the following dimcnsionlcss variables: 

2 

2 \H ^ 

U=K-^aH^^, V = n^cj)H^^^^, £=i^- (H) 
da 

In terms of these variables, equations lf7|)- (|l()|l take the following form (sec the 
Appendix for details) 

= - + X + y + w + m - iv, (12) 

Q = 2e + -i + X - y + W + V - W + {e + 2){V - U), (13) 
= X' + 2X{e + 'i)+y' + 2ye + m{e + l), (14) 
= W + 2W(e + 3)-3V(e + l), (15) 
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where the prime denotes a derivative with respect to iV, where N — hi(a(t)/ao) = 
J Hdt is the number of e-folds. We also note that a/H — a' and = 0'. 

In our discussion we normally assume that the scale factor of the universe before 
inflation is aq, so initially, N = ln(a(t)/ao) = 0. At late times, one has iV ^ 0, since 
a(t) ^ flQ. However, this assumption may just be reversed, especially, if one wishes 
to apply the model to study a late-time cosmology, where one normalizes the scale 
factor such that its present value is oo, which then implies that < in the past, 
since a(t) < ao. 

In a known example of heterotic string compactification ^H^O], one defines, to 
leading order in string loop expansions, A(0) = and ^(ct) = In 2 — ^ e'^ + ct + 

4E^^iln(l - 0-2"^^°'); the latter implies that d^/dcr sgn((T) ^ sinh(cr) < and 

hence Z// > (< 0) for aH < (> 0). For other choice of A((/)), ^(cr) can be different. 
For this reason, as well as for generality, we would prefer to keep both A(0) and ^(cr) 
arbitrary. At any rate, the scalar GB couplings of the above forms may be realized as 
special limits of a more general class of solutions discussed in this paper. 

First we note that the system of equations (|12|I - H15|) may be expressed broadly 
into the following two categories (branches): 

y ^'i + £-]^' -\(J:> + e)U + i(l + e)>V, 
X = -e+\u' -\{\- e)U + i(2 - e)W, 

W = c-^^^-^°\ £ = e(N), U^U{^), (16) 
for V = |yV, where Ni^ is an integration constant, and 

„ , (>V' + 6V)(^-V-2) 
y + 2(2W-3V) ' 

_^_ U'-V {W + ^W-&V){U-V-2) ^ 
2 2(2>V-3V) 
-W-6W + 3V 

£ = , (17) 

for V 7^ |yV. In the first case, for A^ ^ A^o, the time variation of (f) is exponentially 
suppressed, (/)' ~ 0. 

Below we present a number of interesting cosmological solutions (as well as 
construct inflationary potentials), by considering some a prior information about the 
evolution of one of the scalar fields and/or the scale factor. Our discussions also 
generalize past investigations in the literature where comparison is possible. 



^■l. Deleting a -dependent terms 

If we drop the cr-dependent terms, then we find 

C?!)'^ - e^^^'i/^A - = 0, Ck^(/)'0" + - (5 + £)(/)' V3£ + 3 = (18) 

2 

as the set of two independent equations. It is not difficult to see that 0' which decreases 
faster than 1/a^ can lead to an accelerated expansion, that is, e > —1. Note that in 
general A^ depends on the field cj), and its first time derivative, (/> (or cj)'). However, if 
the slow roll-condition ((j)' ~ X^H'^{1 + e) is satisfied at a particular epoch, then in 
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that epoch N primarily depends on the field value. The simplest ansatz is </) = <t)QN+ 
const. Then we find 

=iJoe-P(*/"^''), A(0) = Aoe^P^'^/'*") + Ai, (19) 

where p = (6 + 5Ck^0q)/(6 + Ck^0o)- ^'^^ ^ canonical (f> {C > 0), this solution does not 
give a cosmic acceleration. 

For instance, in an inflationary model motivated by brane worlds, one may 
consider a modified ansatz: AA^ = (l / a) hi{(j) / (/jq) = (l/a)ln0; here 4>o may 
characterize the minimum separation between two moving branes. Hence 

Hi^)^Hj-^^^{e + ,}y'-'^\ A.,.^^, (20) 

where /i = ^k^^qQ^. The universe accelerates for a < —4. Of course, one may 
solve the system of equations (I18II by making an even more complicated ansatz, such 
as (j) ~ 00 tanh(aAAf), which may motivated, for example, by some domain wall 
solutions. The two examples that we discussed above are of simpler kind. 

4-2. Deleting (p- dependent terms 

If we drop all (/)-dependent terms or simply consider a constant dilaton background, 
then we find 

X + y + W U' -{l-e)U -2{X + e)^Q. (21) 

First, as is often the case in most inflationary models, satisfying slow-roll conditions, 
we shall assume that N = cr/cro+ const. Then there exist two classes of solutions: the 
first class of solutions is characterized by 

C(a)=e„e^'-/'^o+ei, V(a) = l-^H(a)2 = ^ii_^^V„e-^'^/'^o, (22) 

where (3 = 1 + 3(5 and 6 = 7K^cr^/2. These are nothing but the GB coupling and 
the scalar potential considered, for instance, in |25| : such terms may be generated 
also using the formalism discussed in The second class of solution may be 

characterized by the potential 

V{a) = (3 - in^H^i' - S)H^, (23) 

where the Hubble parameter is a solution of the differential equation: 

H 2-3^2^2^ • ^ ' 

Especially, for the power-law expansion, H = Hoe"^"^, or equivalently a{t) oc 
{cot + ti)P, we find 

e(cr) = ^oc'^P'^"/""^ + c(l+3P)('^/'^») -I- ^2, (25) 

where = p(l — ^p)/[(l + p)'^^Hq], and and ^2 are arbitrary constants. Unlike 
result H22|l . the GB coupling is now a sum of two exponential terms. 

In fact, alternatively, one can solve the equations H21|) by making reasonable 
ansdtzes for two of the variables; in particular, for X = xq e2"i'^ and U = uq e"^'^, we 
find V{a) = 3mp;i?2(i _ Xqcf^ ~~ Aicr"^/"^), where Ag and Ai take the sign of xq and 

uo, respectively. Similarly, with ai ~ 0, we find V{a) = m^H'^iVo + Vi b-^^VtV^'^). 
The cosmological applications of these potentials will be discussed elsewhere. 
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Figure 1. A non-constant dilaton case. The parameters are chosen as yo = 3, 
(TO = 3/5 (cf equation I3UI ^ and oc e"''^^/-^. There exists a solution crossing 
the barrier w = —1, only if /3 > 5. For /3 < 5, the universe passes smoothly from a 
decelerating phase («) > — 1/3) to an accelerating phase (w < —1/3), but it does 
not cross the barrier w = —1. 



4-3. Late-time cosmology: two scalar case 

Consider that the field a is rolhng slowly such that K\a'\ = \/2(To/7 and kH~^ = 
yJy^V'^l'^ , where (Tq is a (small) constant and i/o > 0. These are good approximations 
for late-time cosmology. Then from H16(l we find 



(20, 

where = 3 + | ^ ± (dg + and 

H oc e('"-^)N A+'"/^ (27) 

where m = 2(yp — ^a^jio^ + 2/0^" '^)- ^ cosmic acceleration occurs for to > 0, which 
is equivalent to the condition V(a) > 2K{a) in the standard model, i.e., without the 
GB coupling. A phantom- like equation of state, namely w < — 1 (or e > 0) may be 
obtained by satisfying C,k?(I)' < 2(j/o — Sctq — 3), or to > 1 if is constant. 

One notes that 4>' = — f + c, where c is an arbitrary constant. The scalar potential 
may be given by 

V{a) e'"'^'" (Ao + Ai e" «fta^'"/4 ^ ^28) 
where /3i = {m — l)y^^^ and /32 = Equivalently, 

V{a) ^ V{a, cj,) = Ao e^^^'' [k' (0^ _ 4c0)] , (29) 

up to a shift in (j). A scalar potential exponential in a (and multiplied by a non-trivial 
function of the dilaton field (j>) may be motivated in string theory, via a study of 
gaugino condensation; in this context, a may be viewed as one of the Kahler moduli 
(see, e.g. EiE)). 
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We can manipulate the system of equations (|17|l so as to obtain the expression 
for e: 

fTo+yo + 3 + W ' 

where W = This class of solution is available only if 7^ const, more precisely, 

if W + 4W 7^ 0. A solution satisfying W + AW < 2(yo - 2(To), or with /5 > 4, is 
inflating for all times (cf figure . Let us also restrict the solution by making the 
ansatz kc/)' = y/2/(e~^^^^^, where [3 is arbitrary. In this case, the EoS parameter w 
is less than —1 when 

n ^ 2 ,/2 ^ 2 yo - 5cro - 3 



If (j) is rolling with a constant velocity, namely k4>' = \/20o/C, then a cx e v?/^^ (<A/mpi) • 
the scale factor grows exponentially with (j). For the power-law expansion, a(t) oc t^, 

we get K(/) = y^plni and Ka = ^J^plnt. 
4.4- Absence of potential, V{a) = 

In this subsection, we study the case where there is no potential for the field a, and 
thus 3^ = 0. This corresponds to the case examined previously by Antoniadis, Rizos 
and Tamvakis [2n| (see also [211123 )i who studied the system (|7|)- (tTnt numerically. 

First, as a special case, we shall consider the solution e = Eq, or equivalently 
a{t) oc (cot + ti)P. The set of three independent equations of motion may be given by 



2 



ci - / e 



2 „2 , _,2_, , W' + (5 + eo)W 



K a i^aH +-^o- =l^ — — ^ , (32) 

6 3(1 + £0) 



.2^/^ „2_ W' + 2(3 + eo)>V 



3(1 + eo) 



where W = Eq = — 1/p and ci is an integration constant. In particular, for a 

logarithmic dilaton, <j>{t) oc ln(cot + ti), and thus (p/H = const = Q, we get 
:i«;V'=cie(i/P"5)N_i^2Q2_ 3(p-^^ 

2 2 ^ 5p - 1 ^ ' 

Let us assume that ci is positive. There can exist some real differences between the 
cases p > 1 and p < 1, and also between p > and p < Q. For 1/5 < p < 1, the 
modulus field a normally behaves as a canonical scalar, even if C > 0, since kQ < 1. 
While, for p > 1, and also for p < 0, since the last term in (|33|l is negative, cr may 
behave as a phantom field, especially for N >0. 

In the following discussion, we analyze a simple and physically motivated case 
where the number of e-folds is proportional to a shift in field cr, namely na = 
^J2aQp) N-V const %. In this case, the set of three independent equations of motion 

X In the V{cf) = case, the field equations are symmetric in — ► cr, except that 5, a- — > —A ^ (or 
U -V) and 7 ^ C. 
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Figure 2. A non-constant dilaton case (cf equation 1341 ') witii «(/!>' oc 
/3(N-No)/2 -poT P > 5, tliere exists a crossing between tlie w < —1 and w > — 1 
phases, when N ~ Nq. The universe is not accelerating for any value of /3 when 
^ ^ ^0- Here wc have chosen ao = 1, but the behaviour of the solution is not 
changed much in a wide range of (tq (> 0). 



reads as 

(1 + e)(2cro + 4(70 = 0, 

where U = K^a'H^^^cr and V = K^cj)' X,^!!"^ . It is not difficult to see that U ~ 
imphes only a decelerating universe; the combined effect of the (scalar) fields is such 
that the EoS parameter w = —1 — 2e/3 = 1, that is, the EoS for a stiff matter. 
However, with U 0, there exist other possibilities; for example, for \U\ > 2ao, the 
universe is accelerating iiU < 0. This acceleration can be super-luminal {w < —1), for 
> 14 > — 4(To; this is consistent with the numerical results discussed in j^O]. Below 
we shall study two more special cases. 

(i) Let us make the ansatz k,{4)/H) oc a^^'^, more specifically, kc/)' = 
e~''('^~^i)/2, where /3 is an arbitrary constant. A solution of this type, with 

/3 > 5, may be used for explaining a late-time acceleration of the universe. More 
specifically, for /3 = 5, we find w ~ — 1 when N < Nq but w > —1/3 when TV > Nq. 
For /3 > 5, there exists a narrow range of A'' where the EoS parameter w falls between 
-1 and -1/3 (cf figure 0). 

(ii) Next, let us consider the asymptotic solution K(f>{t) = ln(a(t))+ const. 
Then we find 



0-0 + 00 ' V 2 V 3 o-Q + 0( 

(35) 



2 3(00 + fTo) 
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Acceleration is possible only if —3 < (70 + 00 < 0, which requires that either 0o < 
or (Tg < (or both). This is essentially an example of phantom cosmology. 



4-5. Small kinetic terms 

Consider the case where the kinetic term for a is absent or it is negligibly small, 
(j/i? ~ §. Then there exist two classes of solutions: the first one is given by 

1 



3 



tt2 

V{a)=c^—e 



2N 



£2 g_2N 

3 



^ e-2N. 



,2N 



(36) 



This is the critical solution with zero acceleration, a (x t, or H = HqC ^, implying 
that V{a) = const. The second solution may be obtained by solving the equations: 



C3 



V = 



w y' + {e-i)y 

3 3(1 + e) 



y + 2ye 



(37) 



In the case the kinetic term (f> is also vanishingly small, 
yields 



3(1 + e) ■ 

i 0, so that W ~ 0, (EH) 



V{a) = TUpiH'^ exp 



(1 - e)dN 



C3 



3 y"(l + e) exp - Ji'^- 



dN 



(38) 



"■PI 



V{a) 



^2 3i74 ' 

A prior knowledge about the time-dependence of the scale factor helps for retrieving 
the functional form of the potential V{a) (or vice versa). For the power-law expansion 
H cx e^"^, we find 

V{ct) = Vo e2^"N + Vi e(i+^°)^, ^{a) = ^eoN ^ g(i-3eo)N^ (39) 

where Vq and ^0 may be (uniquely) determined, but Vi and ^1 are arbitrary. 
Alternatively, 

V{<7) = Vo a-'/" + Vi a-(i/"-i)/"^ (40) 

where a = — 1/eo- This reveals an interesting possibility that during the matter 
domination era, a cx t^/^, the scalar potential may decrease (slowly) as y cx a"^/^ 
rather than the usual result V cx . Moreover, in the presence of matter fields, this 
behaviour can change; in particular, with ci = 0, one finds V cx a"'^'^^"'^ where w 
characterizes the EoS parameter of a barotropic fiuid. 

§ This then represents the case of an almost stabiUzed volume modulus; the product term is 
not essentially vanishing. 
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4-6. Absence of modulus- GB coupling 



Let us consider that U ^ Q (or ^(cr) = 0) and also assume that a is roUing negUgibly 
slow, (T ~ 0. In the absence of scalar potential, V{cr) = 0, we find 
2 



■ exp 



/ 2 



(5 + e)diV 



6 



ci - 3 / (1 + e) exp 



{5 + e)dN 



dN 



(41) 



where ci is an integration constant. These arc the same set of equations as in H18() . so 
wc would not repeat our analysis for this case. With a non-trivial potential. V{a) ^ 0, 
the system of equations reduces to 



ye 



-5N 



Cl 



(y-3)(2y^y) ■ 
2:^3/2 ^ 



(42) 



C0' - e^'i/^A, 



-I - 2Ao, 



where y = Aoexp[— 2 J edN]. Typically, the rolling of a is proportional to the time 
variation of the scalar potential, but with ct = 0, V{a) = Aq/k^. Especially, for the 
power-law expansion, H = H^e'^^^ (eg 7^ —1), wc find 

3(1 + eo) 



=- cie 



-(5+£o)N 



C2e-^^°^\ (43) 

^ + eo 

where C2 is another integration constant. One may consider the possibility that eg > 0, 
or a(i) cx (^oo ~ Cit)" with a < 0, which implies a phantom like "equation of state", 
w < —\. In the case the numerical coefficients c\ and ci are positive, it is possible to 
attain w < — 1 even with a canonical dilaton (C > 0). 



^.7. Absence of GB terms 

By ignoring the GB couplings, so that U — and V = 0, one obtains perhaps the 
most canonical cosmological model with two scalars. If only a is evolving and is 
constant, then of course we find the standard result: y = 3 + e, X = —e. However, if 
(j) is also evolving, then these relations are modified as 

IW 4>" , , 

y = 3 + e, X = -e-W. e = -5 - -- = -3 - ^. (44) 

The implication of the additional field <j> is clear: while the (functional) form of the 
potential is not changed, the kinetic term is modified. Let us discuss the result (|44f) . in 
some detail. At the onset of inflation, e = — 1, one has (j) oc 1/a^. During inflation, in 
at least some region of flcld space, e ^ and hence oc Clearly, a linear dilaton 

background, (j){t) = a + f3t, which implies 4>' = /3/H and cj)" = —PH'/H^, is not 
a solution to the system of equations (|44|l . while, a logarithmic dilaton background, 
(/> oc Int, is clearly a solution. But this last ansdtz implies, in conjunction with H44|) . 
that aif) oc t^^^ ^ which is clearly a non-accelerating solution. 

If the system of equations 144|) is to support a transition between w < —1 and 
w > — 1 phases, then one needs to find a solution for (j) such that its evolution takes 
the ratio (t>" /(t)' from < —3 to > —3 j|. The solution (f) oc with A > 2 inflates 

II The requirement of 4>" /(j)' = {ip/tp) — e < —2 for a cosmic acceleration is not to be confused 
with the slow roll type condition \(j>\/(j>'^ <^ 1 often used in the literature. 
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forever. In particular, for the power-law expansfon, a oc {c^t + ti)" or H (x e we 
get 

It is possible (but not restricted) that </> and a both vary as ncj) oc y^ln^, 
Kcr oc y^lnt; their combined effect may simply be that l^(cr) oc 1/t^ oc a^^/". 

In a matter dominated universe, this yields V oc 1/a'^ |H5I I12j. Of course, in the 
absence of GB couplings, there is no effect like V{a) oc 1/y/a coming from the second 
term in (gU)). 

To summarize, we have seen that, with the standard approximation: H{a) oc 
V^^^{a), K{a) oc H'^i'j) and (j) = const, the universe accelerates for V > 2K with the 
equation of state parameter w > —1. This result is independent of the strength of the 
modulus-GB coupling. If <j) is non-constant, such as </>' oc e~^^^^^, a smooth crossing 
between w < —1 and w > —1 phases is possible for /3 > 4. If the scalar potential, 
V{a), is vanishing, then such a crossing is possible only if is behaving as a phantom 
field and/or the quantity l/l (= K'^a'H'^^^a) changes its sign. Another interesting, as 
well as phenomenologically safe, example is a ~ 0, that is the case with an almost 
constant modulus field. It is possible to attain w < —1 in this case by suitably 
choosing the integration constants (cf H43|) ). even if (p is canonical (C > 0). In the case 
when there are no GB-type corrections, the solutions with a non-constant dilaton are 
particularly interesting. The solution = (poo + 4>i ^"^j with A > 2, inflates forever; 
since w = l + |^, A = 3 corresponds to the case of a pure cosmological constant 
term, wa = — 1- 



5. Second-order system 

Let us introduce the following set of variables: 

7^2 f &Y .^Vicr) 



H " 




(46) 



Note that U and V defined above are different from those in equation (|ll(l ; here they 
do not involve derivatives in their definitions (see the appendix for additional details). 
These provide an alternative way to simplify the study of the model further. Equations 
(f7|l- (|10|) now form a different set of differential equations: 

Q = ?,{U' -V')~Qe{U -V)+ X + W + y -i, (47) 
= V" - U" + {2- e){V - W) - 2{V -U){e' + + 2e) + X + W - y + 2e + 3, (48) 
Q = X' +y' + QX + 2e{X + y) + 3(1 + e){U' - 2eU), (49) 
= W + 2(e + 3)W - 3(e + 1)(V' - 2eV). (50) 

Only equation H48|) is second order in derivative, which can however be discarded since 
any solution following from the other three equations will automatically be satisfied 
by this equation. 
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In a constant dilaton background, for which W = 0, the above system of equations 
simphfy further. There arc broadly three different classes of solutions: 

(z) V = 0, e = e{N), (51) 

(m) V = V(iV), (52) 
(m) V = V{N), e = -1. (53) 

The first branch, (|51|l . corresponding to the case with no dilaton-dependent terms was 
analyzed before in [35]. The second branch, H52|) . implies V cx iJ^ and thus A((/)) = 
const. The first two branches satisfy 

^ + ^-^U' -[e' + e'^ + ^e)U^i + e-y, (54) 

^-^^-^U' -{e' + e'' -e)U = e + X, (55) 
while, the third branch satisfies 

i(V" - U") + 3(V' - U') -^A{V-U)^y-2, 

]^{V" -U") + 2{V -U) = -X +1. (56) 
This is actually the critical case with zero acceleration. 



5.1. Inflating without a scalar potential 

In this subsection we show that in our model it is possible to achieve a cosmic inflation 
without a scalar potential. For simplicity, let us also drop all (/(-dependent terms and 
work momentarily in terms of the original variables (x, u and e) and in the units 
K= 1. Recall that 

u = a^)H^ (57) 

This quantity should normally decrease with the expansion of the universe, so that all 
higher order corrections to Einstein's theory are only sub-leading With V{<t) = 0, 
the equations of motion reduce to the form 

jx^ ^6-6u' + Ueu, e' + e{e + 5) = ^^-^^^^^--. (58) 

2u e 

Fixing the (functional) form of u would alone give a desired evolution for e. In 
particular, for u ~ const = uq 7^ 0, the explicit solution for e is given by 

e = -/3o + P tanh {f3AN) , (59) 

where AA^ = N + C, with C being an integration constant, and 

5w -|- 1 1 / 

Po^^, P^-j2,Ul-2n, + l. (60) 

% Especially, in the case u ~ const, the coupled Gauss-Bonnot term i 5((T)7?.gg = 2>S,{a)H'^(l + e) = 
2iU^H^{l + e) may be varying being proportional to the Einstein-Hilbert term, R/2 = 3H^{2 + e). 
One may need to satisfy the condition \uq\ < 1 at late-times, so that the coupled Gauss-Bonnet term 
is not dominating the dynamics; a phenomenologically safe choice is < 0.2. However, there is no 
such a restriction for early inflation. 
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Figure 3. The EoS parameter w (= —1 — 2e/3) as a function of it^ and A^. For 
My < 0, we get a super-inflationary solution, H /H^ > 0, which corresponds to the 
5 < case first studied in I2UI . Interestingly, for u ~ const, it is possible to get 
a cosmic inflation of an arbitrary magnitude even for 5 > (or u > 0) without 
violating the null energy condition. Of course, ui > — 1 in this case. 



The evolution of a is given by cr' = ± -^(6 + 12eu^/^. The Hubble parameter is given 

by 

H ^Hq e-'^o" cosh(/3AN). (61) 

For Up > 0, both /3p and (3 are positive and /3p > /3, implying that H decreases with 
the number of e-folds, and, in turn w > —\. However, for Mq < 0, we get /? < 0, 
whereas /3q < for —0.2 < < and /3q > for < —0.2. Clearly, with < 0, a 
solution of type H59|) allows the value e > or w < —1 (cf figure (S)). 

5.2. Quintessence/phantom crossing 

Within the best-fit concordance cosmology the present observational data seem to 
slightly favour a dynamical dark energy with Wde crossing —1 from above to below, 
when the red-shift factor z < 0.6 Could it be that our universe crosses the 

phantom divide w = —1, only recently? And could it be due to a GB-type 
modification of general relativity? This is indeed an interesting subject of recent 
discussion on both fronts: observational ISEl and theoretical |3Z| cosmology. Though 
a more detailed analysis of the model, accommodating the effects due to matter fields, 
is needed yet, first indications are that such a possibility cannot be denied. 

First we consider the case where the modulus field a is behaving similarly to the 
dilaton (p, so V ~ U. Then there are again three different classes of solutions: 

{i) X = -e, y = 3 + £, U = 0, (62) 

fji ill ill 

(m) X = -^, y = i+^^ e = ^, (63) 

+ In the presence of matter fields, Q,m > 0, one actually defines an effective EoS for dark energy, Wgjf, 
whose value can be (somewhat) different from that of the dark energy EoS parameter, w = Pde/ Pda- 
Nevertheless, if a (fundamental) theory is to support the crossing of phantom divide, w = —1, then 
a similar effect is expected in the absence of matter since a gravitationally attractive source can only 
make the dark energy source less squishy or the EoS parameter more positive, w^f f > Wde, ■ 
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(m) A- = 1, y ^2, £ = -1. (64) 

The first branch is the standard one, implying that \{4)) = = S,{<^)^ while the third 
branch gives zero acceleration, so it is the second branch which is most interesting. 
For H oc (C(o')) , so that U ~ const, the effect is similar to that of a pure 
cosmological constant term, £ = 0. Instead, iiH (x a~ then the universe 

is accelerating for a < 2, whereas a = 2 marks the boarder between acceleration and 
deceleration. Typically, for the solution U ~ ciN + C2 (with > 0) wc find £ < for 
< —Cijcx and £ > for > ~C2/ci; where = ln(a(t)/ao), with oq being the 
present value of a{t), and 1 + 2 = c'^^. For instance, if C2/C1 — 0.47, then the EoS 
parameter w is less than —1 for z < 0.6. 



5.3. Deceleration/acceleration phase and crossing of phantom divide 

Consider the possibility that the coupling functions ^(cr) and A(0) both vary inversely 
with H^, so that U = uo and V = vq. For uo = vq (which again implies a trivial GB 
coupling), there are three different classes of solutions. The first class of solution is 
given by 

£ = -3 + \/6tanhV6A7V, cx e~ sinh 2\/6 AN, 

^ ^(6wo- l)tanhV6AiV~3(5i;o- 1)1 , (65) 



7(7 



k2 



((j)'^ = --^vo (2%/6tanh Ve AiV - 5 

where A = N + Ni. This solution yields only a standard inflation for which £ < 
(cf figure 0J. The transition from the deceleration phase to the acceleration phase (or 
vice verse) is a natural ingredient of the above solution. In fact, we would not have 
gained such an insight into a functional from of e (and the EoS parameter w), have we 
had dropped the GB couplings from the action, which corresponds to the case vq = 0. 

For |AA^| > 1, since |tanhV6AA^| ~ 1, both the fields a and are rolling with 
almost constant velocities, viz, a = ctoV^A^, where ctq = 'mpi^ (2 — 9wo)/37- The 
scalar potential and the modulus GB coupling take the following form 

V{a) = Vq c~ ^(■"/""'i sinh 2 (c + a/ao) , 

^(a) = eo c^(^/-°) [cosh (c + a/ao)]- ' , (66) 

where c is an arbitrary constant. A similar expression may be given for A((/)), in terms 
of (j). 

The other two classes of solutions may be described by the following system of 
equations: 

X = -e-W, y = 3 + s, (67) 
where, as above, W = and 

2 6^0 6 y V ^^0 / vo Vq 

In the constant dilaton case, (p' ~ 0, wc find = — 1 or 0. Wc assume that we are 
dealing with a canonical dilaton, so W > 0. For W/Vq > 0, the both roots in (|68|l 
give the EoS parameter w > —1. However, for W/Vo < 0, the positive root can give 
w < -1. 
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Figure 4. The universe enters into an accelerating phase (e > —1) when 
AiV ~ 0.22, or H/Hq ~ 0.42 (cf equation ^^). For AA^ > 1, the field </> is almost 
constant, with the deceleration parameter q = —1 — e ~ —0.45 and w ~ —0.63. 




Figure 5. The equation of state parameter w (= —1 — 2e/3) with ft(^' = 
y 2^ e-''(N-No)/2 ipo/vo = 0.2 (cf the negative root in equation 15^1. 

Only a small value of /3 seems to be reasonable, otherwise ui can be imaginary for 
AN < 0. 

To simplify the analysis further, let us make the ansatz Kcj)' = ^ g-/3AN/2^ 

where /? is arbitrary. Then, for (f>a/va > 0, in the regime AA^ ^ 0, such that (p' /4>q — » 0, 
the positive root in H68() yields w ^ —I (except for /3 ~ 0), while the negative root 
yields w ~ —1/3 (see figures |S1 and O . On the other hand, for W/uo < 0, the 
positive root supports a phantom phase (w < —1), while the negative root solution 
is not accelerating, as it approaches w = —1/3 from above (see figures [3 and (HJ . It 
is generally for uq < (or equivalently < or 5 < 0) that one of the branches 
can give a non-standard (w < —1) or super-inflationary phase satisfying p + p < 0, in 
consistent with the numerical results in [201 HH] • 

For the power-law expansion, a{t) cx t", so that H = a/t, and (/)' = const (or 
cf> (X \nt+ const), we require (j>o/vo < 0.30 or > 29.69 so as to keep e real. We shall 
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Figure 6. The equation of state parameter w as in figure|S] but for the positive 
root in equation I68i . For /3 > 0, the universe passes from a decelerating phase 
to an accelerating phase, but docs not cross the phantom divide, w = —1. 




Figure 7. The equation of state parameter w with mj)' = y No)/2 
and 4>o/'t-'o = —0.2 < (cf the negative root in equation 1(381 1. Here w may start 
from a value (well) below —1 but it never crosses it, leading always to a phantom 
phase. 



be interested in the case where vq is small, so that the GB contribution is only sub- 
leading. For instance, if W ~ —O.Ivq > 0, then e ~ —1.06 or 0.09; only the second 
branch is accelerating which gives w ~ —1.06. 

In the Mo =/= Vq case, so that the GB coupling is non- vanishing, the first branch, 
H65|l . is modified as 

7^^' ' = 4 [(7wo - wo + 1)(3 + e) - 3(6uo - ^'o - 1)] , 

Hi 

V ^ 7n%i{vo - uo + l){3 + e) (69) 

without changing the forms of W and e. As for the second branch, (|68|) , the expressions 
of 7(t' and V{a) are now more complicated and not illuminating at least to write down 
here. However, a simplification occurs in the case that H = Hq e^°^ or equivalently, 
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Figure 8. The equation of state parameter w as in figurc|7| but for the positive 
root in equation I68i . For any value of f}, the universe is not accelerating, since 
in this case w always takes a value larger than —1/3. 



a cx (cot + ti)" (a > 0) or a oc (to 
the scalar fields are given by 



t) " (a < 0). Especially, the time variations of 



/2 ^ 



CiC 



-2(3+£o)N 



cie 



3(1 + £o)£oUo 

3 + eo 

-2(3+eo)N _ 3(1 + £o)goUo 



+ £o(l - £o)(uo - Vo) - £o 



3 + .„ ■ <™' 

Here either value of £o is possible; for £o > 0, (/> may behave as a phantom field, 
especially, if > and/or ci < 0, unless that uq (or ^(ct)) is negative. 



5.4. Late-time approximation: two scalar case 

Finally, as a good approximation at late times, consider that one of the fields, mainly 
(7, is rolling very slow such that h\cf'\ = ^J2xq/^ and K,H~^{a) = V^~^/^(cr), where 
Xq is a small constant, while j/o > 0. In the case of a vanishing effective GB coupling, 
U = V, the system of equations reduce to 



C0" 



«2e(a) 



ci e 



^2(3-yo)N _ yo(xo + Yo - 3) 

3(yo-2)(3-yo) 
(3 - xo - yo) , £ = yo - 3. 



(71) 



The universe accelerates for yo > 2. An especially interesting range for yo is 
2 < yo < 3, for which the universe naturally undergoes a cosmic acceleration and 
the couplings ^(tr) and A((/)) both decrease exponentially with the number of e-folds, 
N. 

In the general case, there are five different branches of solutions: except the first 
branch, which is giving by 

e(a)=0, e '"-^ 



C(p = — 



ci e 



xq + yo 

.(2+£)N ^ (2xo-yo)(3-xo-yo) 
2xo + 5yo 



(72) 



On compatibility of string effective action with an accelerating universe 



22 



1 



-6xoN/(xo+yo) 



C2 



^6xoN/(xo+yo) 



xo - yo 



C 2 

2^' 



dN 



k2A(0) 

their expressions are lengthy and not iUuminating at least to write down here. 
However, especially, for e = epj they merge to give a single solution 



m 



yo(l - £o) - 3(1 + Eo) 



Xo — 3Ai(5 + 3eo) 



, ^0 + Ci c" 



-2£nN\ 



(73) 



1 



(Ao + Aie-(5+3-«)N^;^^g-2soNj 



where 



Aq and Ai are integration constants, and 
3a;o + (a;o + yo)eo 



3£o(l + eo) 
_ (3 + Eo) [(3 + a;o + yo)eo + 3 + Sxq - yo] 



(74) 



Kcr so as to retrieve a functional form 



3eo(5 + £o)(l + Eo) 

with £o 7^ ^1- One may replace N by 

for ^(cr). Especially, the cases Eq = 0,-5 (and also Eq = ~5/3) need to be treated 
separately, for which new terms appear with iV-dcpendcnce. For instance, for eo = 
(i.e. H = Hq), we find 

'yo - 3 



Xo + ci e 



- 5N 



A(0) 



1 



'-Hr 



An e 

15 

Ko - 2a-o7V] 



1 



3 - yo)N 



(75) 



With no surprise, this solution behaves all the times as a cosmological constant term. 



-1, since H = Hq and V{u) 



"-PI 



for which w 

We conclude this section by providing an overview of the different solutions 
discussed above. In subsection 5.1 we showed the existence of an inflationary solution 
with a vanishing potential and for the modulus-GB coupling ^(cr) oc . This result 
is quite interesting. It is also obvious from the above analysis that solutions with 
trivial scalar-GB couplings {U = = V) are not equivalent to solutions where these 
coupling are (negligibly) small, at late times. The simplest way to visualize this is 
to set only U = V, i.e., the dilaton </> is behaving as the modulus a. Unlike in the 
U = = V case, we now find three different kinds of solutions, one of which, namely 

can support a crossing between w < —1 and w > —1 cosmologies. When A(0) and 
^(cr) both grow inversely with the squared of the Hubble parameter, iJ^, then, for the 
first branch, the EoS parameter w can take a value less than —1/3, namely w ~ —0.63 
at late times. Either both of the other two branches support an inflationary phase 
with w > —1, or only one of the branches supports an inflationary phase, but with 
w < —\. In the latter case, the second branch is not accelerating since e < — 1. 

As in the standard picture, with H cx V'^/'^ and a ex H(a), the universe accelerates 
for = K^-'^jT^ > 2, even \iU — V ^Q. This result is independent of the ratio &/H. 
It is the case for < 3, the function U decreases with the expansion of the universe. 
In the general case, for which ^ V 7^ 0, there are five branches of solutions, which 
simplify, for instance, if we take the limit where the slow-roll type variable e is constant. 
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6. Runway dilaton and matter coupling 



In a more realistic scenario, V{a) in ((^J may be replaced by a general form of the 
potential, namely 

V{a)-*V{4>,a) = \{<t>)V{a), 

which is now a function of both the modulus a and the dilaton (p. In our analysis 
we considered only the simplest case where </> was expected to evolve (relatively fast) 
towards its finite value 4>m- In this limit, \{4>) ~ const. 

Generally, (p (the spin-0 partner of spin-2 graviton) is coupled to Einstein term and 
thus, while going from string frame to Einstein frame, it may also couple to matter 
fields. In the heterotic string theory studied by Antoniadis, Gava and Narain 
there is no coupling between the field a and the Ricci-scalar term, so a couples 
to gravity only via spacetime curvature tensors or Riemann curvature invariants. 
Thus a conformal transformation does not necessarily involve a cr-dependent term, 
but only a dilaton dependent term. In addition, a generally has a non-zero vacuum 
expectation value and hence a non-zero mass |2()| : it is thus conceivable to expect 
only a minimal coupling between a and the matter. Furthermore, in a de Sitter 
(accelerating) universe, sypersymmetry has to be broken in general. If so, dilaton 
acquires a large enough mass (m^ > 10~^eV) after supersymmetry breaking. Within 
this scenario, 4> takes an almost constant value and hence it decouples from matter, 
so called least coupling principle |38l I39j . In any case, let us assume that is coupled 
to matter fields non-minimally: 



- y (a) - ^ (V0)2 + 1 TZIb + AH^) [PM + PR) 



.(76) 



where / = [A(0) — S^{a)]. The equations of motion that describe gravity, the scalars, 
and the matter and radiation are given by (in units k = 1) 



H^ = l + ^^' + V- -SH-'f + A^ipM + pr) 



a 
a 



7ct' + C4'''-V + A'pr + ^A^pM + ^ i?' [hJ + f + 2fH/H 



1 



C((/) + 3i?0) = - A, ^A' PM + ^ A, ^ U'cB 



(77) 
(78) 

(79) 
(80) 



In the standard scenario of no matter-scalar coupling, e.g. the model of quintessence, 
A{(j)) = 1. For ^4(0) 7^ 1, the fluid equation of motion for matter (or radiation) is 



dpi , 1 dp, 
oa a 0(p 



(81) 

where i = A/(matter) or i?(radiation) and a = The implicit assumption is that 

matter couples to A^{(f))gfj^^ (with scale factor a) rather than the Einstein metric g^i, 
alone, and pR oc where a = aA{(j)) so only enters the equation of motion. 

We can practically drop the pr term, at late times, since most of the radiation energy 
in the universe is in the cosmic microwave background, which makes up a fraction of 
roughly 5 x 10^^ of the total density of the universe. 

In the case A{<f)) = 1, one can solve the system of equations (I77ll - (|81|) . for example, 
by making an ansatz for the evolution of one of the scalars and/or the scale factor of 
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the universe. Note that in this ease the energy eonservation equation (|81|l is easily 
expressed as a differential equation, namely il- + 2£f7j + 3(1 + Wi)^i ~ 0, where 
f2j = pi/3H'^, Wi = Pi/ Pi and Q! = dil/dln[a(t)]. In the ease ^ 1, the extra 

(matter) degree of freedom added to the model complicates the problem, in which 
case the field equations may be solved numerically |40| . 

Especially in the case where (p evolves considerably slow with time, A{(j)) is not 
essentially a constant, even at late times. In this case, a question may be raised as to 
whether the equivalence principle (universality of free fall in a given gravitational 
field) is violated due to a runway dilaton j^Hl EHli leading to time variations of 
some of the fundamental constants, like the fine structure constant ae = e^ /he. 
For metrically coupled matter theories, one may also raise a question regarding the 
time variation of Newton's constants under post-Newtonian approximation |2J . But 
the constraints on gravitational actions of the type considered here, under (post- 
)Newtonian approximations, may be less stringent as they only require [421 143j (in 
units K = V) 



&H^-^ < 10"^ pH^—^ < 10'\ (82) 
do" d<p 

In our notations, introduced in equation Hll|l . these read UH < 10^^ and VH < 10^^. 
This result does not rule out any scalar- GB couplings. Within our model, both U and 

V may be exponentially small (close to zero). In any case, ifU ^ H^^^^j^ ^ C'(l) ^^^d 

V ~ -ff^^l^ ^ C(l), then we would require that fi, < 10"^ (in units k = 1). 

If X{(j)) and ^((t) are nearly constants, especially, after inflation, which we assume 
to be the case here, then action H76(l may be compared to the gravitational action 
considered in [SH| (cf equation (2.4)), under the replacement, V{(t) — > V{(j),(T) = 
X{(l))V{a). For most of the solutions that we have presented in this paper, after a 
certain number of e-folds of expansion, the rolling of cj) could be small as it is given 

by (j)' = ^J^^H ~ e~^^ with (3N > 0. In a specific model considered in |38| . 

X{(f>) = Aoo(l + bxc~'^'^), for ^ = (j)' < 0.7, the level of variation of fine structure 
constant, \dlne^/dt\ < 10~^^yr~-'^ is much smaller than the current best limit on the 
time variation of e^, namely \d\ne'^/dt\ < 10~^^yr~^; thus, a small rolling of is 
physically viable and phenomenologically acceptable. 

One may be more concerned here with a time- variation of the (common) modulus 
a. We take rr = as a conservative case, for which V{a) is also constant (cf 
equation (|42|l ). Because of phenomenological reasons it is desirable that the rolling 
of a is small, especially after inflation, which then allows only a slow rolling of the 
potential. Instead of the condition ct = 0, we expect k& < H to hold at all times, 
i.e. ^2X < y/^. This last condition is satisfled by all our solutions discussed in the 
previous sections. 

In has been known that even if the tests of general relativity in the solar 
system (and binary pulsars) set strong constraints on scalar-tensor models of the 
type considered in this paper, they may differ significantly from GR at high red- 
shift, z > 0.7 |44[ I42| . possibly affecting the standard scenario of structure formation. 
In turn, the scalar-GB couplings, in addition to the matter-dilaton coupling, may 
be constrained only by taking into account cosmological and solar-system data 
together 02]). Such results are not known yet; see also the recent paper which 
discusses on cosmological and astrophysical constraints of Gauss-Bonnet dark energy. 
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For A{(t)) ~ 1 and / ~ 0, one has </> oc (cf equation (|80|l 'l: with the large 

scale factor, the rolling of 4) can be small. In fact, the universe has merely doubled in 
its size between radiation domination era {t ^ 10^^ sec) and today {t ~ 10^^ sec) or 
A'^ = 1 + z ~ 2, where z is the red-shift factor. So, unless we do not demand a large 
number of e- folds (between radiation domination and present eras), we find no reason 

to expect a large shift in 4>, since A0 = J \/WHdt ~ ^J'^^ AN . It may well 

be that the field (j) has evolved by a small amount; provided A,^/A < 10~^, or more 
mildly < 10~^ jl^l, for today's value of (f), any possible violation of equivalence 
principle (or deviations from the standard general relativity) may be small enough to 
easily satisfy all cosmological constraints. 



7. Conclusions 

In this paper, we explored the cosmological implications of the modulus-dilaton 
dependent one-loop corrections to the gravitational coupling of the heterotic string 
effective action, by presenting various exact cosmological solutions. These corrections 
comprise a GB invariant multiplied by universal non-trivial function of the common 
modulus field a and the dilaton field 0, and thus they may have an important role 
in explaining inflation, as well as the recently discovered acceleration of the cosmic 
expansion. The solutions represent a flat four-dimensional FRW spacetime where the 
Hubble expansion parameter is a smoothly decreasing function of the number of e- 
folds, N, while remaining non-zero at all times. As such, they exhibit the type of 
branch-changing solution, i.e., deceleration/acceleration regime, that is a prerequisite 
for the successful incorporation of the early inflation as well as the currently observed 
deceleration/acceleration phase in the expansion of the (late) universe. 

The main purpose of this paper was to point out the importance of a'-corrections 
for explaining the cosmic acceleration problem attributed to "dark energy" or 
gravitational vacuum energy. It may well be that the evolution of a cosmological 
(FRW) background from the string perturbativc vacuum leads to a regime of cosmic 
inflation in which the strengths of the coupling constants associated with the string 
loop effects can grow with time, though very slowly, while the higher curvature 
contributions to the effective action are more suppressed than the Ricci scalar. When 
both the common modulus filed (a) and the dilaton (</>) take almost constant values 
at late-times, the model reduces to a standard picture. 

We demonstrated that, both for slowly running dilaton phase and constant dilaton 
phase of the standard scenario, the model can lead to a small deviation from the 
u' = — 1 prediction of non-evolving dark energy; some of the branches, specifically 
when the parameter S in the one-loop effective action is negative (or ^(cr) < 0) can 
support a non-standard (phantom) cosmology (w < —1), while the other branches 
support a standard cosmology with no-acceleration {w = —1/3) or acceleration 
(— 1 < w < —1/3) or deceleration {w > —1/3). 

In our analysis, we have assumed that there is no potential for the dilaton in the 
loop-corrected perturbativc string effective action. However, it is quite possible that 
for the formulation of a complete and realistic string cosmology scenario, one needs 
to take into account the quantum back reactions of loops and radiation as well as an 
appropriate non-perturbative dilaton potential. 

We assumed that in our model the contributions arising from two (and higher) 
loop terms are suppressed as compared to those from one-loop terms. In this sense, 
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the result discussed above may persist in the full theory under certain assumptions 
about the moduli dependence of loop corrections to higher derivative terms. A 
key point is that a field dependent (and thus space-time dependent) GB coupling, 
which is completely reasonable in the four-dimensional effective action, can make the 
observability of various cosmological parameters, including the Hubble expansion rate 
of the universe, more achievable. 

A more detailed analysis of the model, mainly by introducing matter fields, is 
needed yet, first indications are that a GB-type modification of gravity is a healthy 
proposal for explaining the cosmic acceleration problem. 
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Appendix: Field equations as a system of differential equations 

In order to express equations l(7)l-l(TU)l as a set of differential equations, we define the 
following variables: 



V 



One also notes that / = A(^) — ^{a), and 

; • •• -nd^A dA ■■ 

^='^V+d^'^' (^-^^ 

and a similar expression for ^. A simple calculation yields 



A2 



Q'^i - 2ew a^A ^2 - 2ew 

dV I _y' + 2ey 
da X ' 



(A.3) 



where the prime denotes a derivative with respect to N , where N = J Hdt = 
\n{a{t)/ao). Note that ^ = = H The equations of motion, equations Q- 

P0|l . then take the form 



= -3 + y [-fx^ + Cz^j + n^y + 'in\xu - zv), (A.4) 



(A.IO) 
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= 2e + 3 + ~ + K^{£ + 2){zv - xu) + — (73;^ + Qz^) ~ n'^V, (A. 5) 

= — — (7x2) ^ ^2(^ ^ 3)^^2 _^ ^2^y, _^ 2ey) + 3k2(1 + e)xu, (A.6) 

= y ^ (C^^) + + 3)7z^ - 3a^^(1 + (A.7) 

Finally, by defining the new variables 

^2 ^2 
^" = -^7x2, y = n'^y^ y^=^^x^^ U = ti^xu, V = k'^zv, (A.8) 

we arrive at equations (fT^ - (fT5|l . 

Alternatively, equations (|7j)- (fTn|l may be expressed as a set of first- and second- 
order differential equations. To this end, one modifies the definitions for u and v in 
equation (|A.ip . namely 

u^^{a)H^, u = A(0)ij2. (A.9) 

A simple calculation yields 

d£ u' - 2eu , n u" - 2e'u - 2eu' , , „ , u' - 2eM 

"5-^" = ' = 5 (x +2e.T) , 

aa X O'^a x^ x"^ 

dX . v'-2ev ^^A „2 v"-2e'v-2ev' ,, „ ,w'-2ew 

ocp z O'^cj) z^ z-^ 

By substituting these expressions, along with the first, second and fifth relations 

defined in (|A.3|I . back to equations l(7|l- pU|) . and again defining new variables as in 

IjA.Sp . except that U = k^u and V = k^v, we arrive at equations l|T7|) - (|5n|) . 

Finally, we find several fixed points for the system 1|A.4(I - (|A.7|I . characterized by 
x' = y' = z' ^ 0, and u' = or v' = 0. The fixed-point solutions are specific to the 
choice of the potential or the couplings ^(cr), X{(f>), other than the coupling constants 
7 and C; below we will make a canonical choice, that is, C = 2j/3. 

For the system of equations ljA.4p - (|A.7p . there exists a (global) fixed-point, as 
given by (x, y, z) = (0, 3, 0), which corresponds to a stable de Sitter phase with e = 0. 
In particular, in the case of a vanishing potential, so y = 0, the fixed-point solutions 
are given by 

• (A) (x, z, u) = (xo, 0, 1/xo — 7x0/6), with an arbitrary v, and e = — (6 + 57Xq)/(6 -|- 
7Xq), which gives a (stable) de Sitter solution only if 7 < 0. 

• (B) ix,z,u) = ( -3"o±V67+9« |^Q^^^^^ ^j^j^ arbitrary w and e = -3 - ^^=f , 
which gives a de Sitter solution for uq < —^^2^/3. 

• (C) (x, z, v) = (0, zq, ^ — ), with an arbitrary u and e — — Vit''~Jl ■ This branch 
gives an accelerating solution only if 7 < 0. 

• (D) (x, z, v) (0, ^''"^^■^^"'^^ -,vo), with an arbitrary u and e = -5^o . 
Depending upon the sign of vq, one of the solutions gives a stable de Sitter phase. 

If J/ 7^ 0, then, for the purpose of finding fixed point solutions, one needs to 
specify the functional form of V{(j) and/or the coupling ^(cr) or X{(j)). As a simple 
case, consider V — Voe~'^°', which implies that y' + (2e -I- Xx)y = 0. The fixed point 
or critical solutions arc then given by 
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where z = zq e~^^ and u = u{N). 
.(F)a-=f, y^^, z = 0, e = -1, « = 



where v = v{N). 



where z = zqB , u = ^ (y " ^ ^ with c being an integration constant, and 

V = v{N). 

If one restricts the couphng f (ct) to be ^(ct) = e^*^, or equivalently u' + 2(1 — 
/i/A)u = 0, then the above critical solutions are further constrained; we now have 

• It = „ ^ 

9+272:^ 

• A = /i 

respectively, for the branches _E, and G. It is obvious that the choice A = /i is 
special. 

For C, ~ 27/3, the system of equations (|A.4p - ljA.7() can be reduced to the following 
two sets of equations: 

dx 1 du 2e + 72;^ 
dN=^udN + ^^+''^'-'^' 

+ 3ux^0, z==0, (A.ll) 



and 



dz _ 6y - 18 + 37X^ - 1672^ + 18ua- _ 9(i; - 72) 
d/V ~ 47Z ~ 27 ' 

(i.T ^ 70:^ z a; dit 772 dz 

dN ~ ~ ^ udN ^ um ^ ~QudN' 



For F = Vo e , from the branch (|A.lip . wc obtain 



(A.12) 



. . viZEM ,,,,747^ ^ 27 

3A 27 ^ °^ 3A ^ ' 

and a; = — 3^/7. In addition, if we choose ^ = ^0 e'*'"', then we require A^ = 27/3 and 
hence du/dN = dx/dN = 0. The stability analysis of the critical solution a oc (t + ii)" 
and a cx \n(t + ii) would be similar to that in |25j . Further implications of the model 
will appear elsewhere. 
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